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(N ■ ABSTRACT 

Giant flares from soft gamma-ray repeaters (SGRs) are one of tlie most violent 
phenomena in neutron stars. Quasi-periodic oscillations (QPOs) with frequencies 
^ ! ranging from 18 to 1840 Hz have been discovered in the tails of giant flares from 

two SGRs, and were ascribed to be seismic vibrations or torsional oscillations of 
magnetars. Here we propose an alternative explanation for the QPOs in terms 
^ ! of standing sausage mode oscillations of flux tubes in the magnetar coronae. We 

show that most of the QPOs observed in SGR giant flares could be well accounted 
for except for those with very high frequencies (625 and 1840 Hz). 
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O ' 1- Introduction 

I> 

^ , Soft gamma-ray repeaters (SGRs) are neutron stars that exhibit sporadic burst activities 

O ■ most prominently in soft gamma-rays (Norris et al. 1991). Besides normal bursts with energy 

E ~ 10^^ ergs, enormously energetic giant flares (with E ~ 10^^ — 10^^ ergs) have also been 

<> '. observed, for example, from SGR 0526-66 in 1979 (Mazets et al. 1979), from SGR 1900+14 

in 1998 (Hurley et al. 1999; Feroci et al. 1999) and from SGR 1806-20 in 2004 (Hurley 
et al. 2005; Terasawa et al. 2005; Palmer et al. 2005). The giant flares generally start 
with an initially rising spike lasting ~ 1 s and then evolve into a decaying phase which lasts 
hundreds of seconds. Theoretically, SGRs are thought to be magnetars, neutron stars with 
surface magnetic fields strengths of ~ 10^^ — 10^^ G (Thompson & Duncan, 1993; Duncan & 
Thompson, 1994; Thompson & Duncan, 1996, 2001), or neutron stars with normal magnetic 
fields accreting from a disk formed from the fallback in supernova explosions (e.g., Alpar 
2001). 

During the three SGR giant flares mentioned above, quasi-periodic oscillations (QPOs) 
were identified both in the initially rising spike and in the decaying tail (Barat et al. 1983; 
Terasawa et al. 2005; Palmer et al. 2005; Israel et al. 2005; Strohmayer & Watts 2005; 
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Watts & Strohmayer 2006; Terasawa et al. 2006). These QPOs seem to provide independent 
evidence for superstrong magnetic fields in SGRs (Vietri, Stella, & Israel 2007). In the 
most popular models, they are explained as global seismic vibration modes of magnetars 
(Hansen & Cioffi 1980; Schumaker & Thorne 1983; McDermott, van Horn & Hansen 1988; 
Strohmayer 1991; Duncan 1998). Other explanations have also been proposed. Levin (2006) 
argued that the QPOs may be driven by the global mode of the MHD fluid core of the 
neutron star and its crust, rather than the mechanical mode of the crust. Following this 
idea, Sotani et al. (2007) recently made two-dimensional numerical simulations, and found 
two families of torsional Alfven oscillations which may explain some of the observed QPOs. 
Coupling neutron star's elastic crust and fluid core, Glampedakis et al. (2006) used a simple 
toy model to provide explanations for some of the QPOs observed. Besides, magnetospheric 
current variation induced by the crust oscillation was proposed for the modulation of X-ray 
flux by Timokhin, Eichler, & Lyubarsky (2007), and they showed that radial oscillation with 
an amplitude of one-hundredth of the neutron star radius can account for the observed QPO 
flux fraction. 

In this paper we propose an alternative explanation for QPOs observed in the SGR giant 
flares. We assume that part of the plasma ejected during the giant flares are trapped by the 
magnetic fields, and then form magnetic flux tube structures, similar as in the solar corona. 
Oscillations of magnetic loops or tubes in the solar corona have been observed by a lot of 
authors (see Aschwanden et al. 1999 for a review). We argue that tube oscillations in SGR 
magnetospheres may give rise to some of the QPOs observed during the giant flares. We 
discuss the plausible tube oscillation modes in the SGR magnetosphere in §2, and compare 
the possible oscillation frequencies with observations based on the flreball model of Thompson 
& Duncan (1995) in §3. The possible implications of our model are discussed in §4. 



2. Magnetic flux tube oscillations 

In solar physics QPOs with periods of several minutes in coronal loops have been de- 
tected and successfully interpreted in terms of MHD waves (Aschwanden et al. 1999; Nakari- 
akov et al. 1999). Following this idea, we suppose that these MHD waves could also be 
excited in the magnetar coronae when starquakes of a magnetar shear its external mag- 
netic fleld, which becomes nonpotential and threaded by an electric current (Beloborotov & 
Thompson 2007). As pointed out by Robets, Edwin & Benz (1984), MHD tube oscillations 
in the sausage mode can modulate the cross section of the tube, and hence its density and 
radiation flux, while oscillations in the kink mode do not change the loop density in the flrst 
order and thus cannot modulate the plasma radiation (Aschwanden et al. 1999). So we rule 
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out the kink mode oscillations for the QPOs observed in SGRs. Furthermore, observations 
of the QPOs show that they can last hundred to thousand cycles of the oscillation periods 
(Israel 2005; Terasawa et al. 2005; Palmer et al. 2005; Strohmayer & Watts 2005; Watts 
& Strohmayer 2006). Accordingly, we exclude the propagating sausage mode oscillations 
because of their rather short damping timescale (about tens of the oscillation period), and 
consider only the standing sausage mode oscillations. 

The flux tube can be roughly described as a cylinder of radius a and length L. The 
magnetic field strength, temperature, mass and number densities are denoted as -Bq, Tq, 
Po, no inside the tube, and Be, Te, pe, ''^e outside the tube, respectively. Magnetoacoustic 
oscillations of a magnetic tube have been studied thoroughly, and the derived oscillation 
periods Tfast for the fast standing sausage mode (Wentzel 1979; Edwin & Roberts 1983; 
Roberts et al. 1984), and Tgiow for the slow standing sausage mode (Edwin & Roberts 1983; 
Roberts et al. 1984) are expressed as follows: 



Tfast 



^ = "^^'^M-^^y^' ^ 2.0 X 10-5^ s, (1) 

Tsiow = — = — [1 + i^n^' ^ 0.015^ s, (2) 

where Ck, Ct, Cg, vx are the kink speed, slow magnetoacoustic speed, sound speed, and 
Alfven speed, respectively (Roberts 2000), Lg = L/10^ cm, a^ = a/10^ cm, n28 = n/10^^ 
cm~^, Bi4 = B/10^'^ G, and Tiq = kT/10 kev, j is the number of nods and j = 1 for the 
fundamental mode. In deriving Eqs. (1) and (2), we have made the approximations that (i) 
the sound speed in the magnetar magnetosphere is much smaller than the Alfven speeds, 
i.e., Cg <^ Va, so the slow magnetoacoustic speed Ct is close to the sound speed cq inside the 
flux tube; (ii) the magnetic field strength inside a magnetic flux tube Bq is comparable to 
that in the environment B^, whereas the plasma density po inside the tube is much higher 
than Pe in its surroundings. Thus, we have Ck — V^fA- We also ignore the effects of 
gravity on the flux tube because of the super-strong magnetic fleld of magnetars. With the 
typical values for the parameters of the magnetar tube, the above equations indicate that 
the oscillation frequencies in the fast standing sausage mode are too hight to be compatible 
with the observed QPOs in SGRs, and we are left with only the slow standing sausage mode 
oscillations. 



3. QPOs in giant flares 

In this section, we investigate how well the slow sausage mode oscillation frequencies 
can match the QPOs frequencies observed in SGR giant flares. The QPOs observed in the 
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initial spike phase and in the flare tail phase are discussed separately. 



3.1. QPOs in the initial flare spike phase (t < 1 s) 

Geotail spacecraft mission was originally aimed to study the structure and dynamics of 
the tail region of the magnetosphere of the Earth. With it 50 Hz and 48 Hz QPOs were 
detected at t = 45 — 175 and 430 — 567 ms after the onset of the giant flare in SGR 1806—20 
in 2004 (Terasawa et al. 2006). The oscillation periods (r ~ 20 ms) are similar to the period 
r ~ 23 ms of the QPOs from SGR 0526-66's giant flare in 1979 (Barat et al. 1983), but no 
similar phenomena have been seen in the SGR 1900+14 flare. 

The physical picture in our flux tube oscillation scenario is as follows: after the onset of 
the giant flare, hot plasma was ejected into the magnetosphere above the surface of the SGR. 
Because of the super-strong field strength of the magnetar, the plasma could only move along 
the field lines, from one to another footpoint of each field line on the star's crust, and could 
not get out of the confined structure of the field. This is why a fireball is formed. We assume 
that the magnetic tubelike structure(s) were subject to various types of oscillations excited 
by the turbulence at the footpoints. For the slow sausage mode oscillation, the period is 
related to the initial e-folding rising time t^se = L/cs, which is 9.4 m^j in the giant flare from 
SGR 1806—20 according to Geotail spacecraft observations (Tanaka et al. 2007), i.e., 

2L 2L 

i"siow = ^ ~ ^ ~ 2trisc ~ 18.8 ms (3) 

JCt JCs 

for j = 1 (we do not need to consider the gravitational redshift effect in TsIow here because 
the observed rising time trise has already included the redshift factor). This period is very 
close to the observed value ~ 20 ms. 

Unfortunately, there was no QPO detected in the initial rising phase of the giant flare 
in SGR 1900+14 with rising time tnse ~ 3.1 ms, and no rising time was measured in the 
1979 March giant flare event with an initial ~ 43 Hz QPO. So in the three giant flares ever 
detected we have only one event on SGR 1806—20 to examine our explanation for the QPOs 
in the initial spike phase, which should be testified by future detailed monitoring of SGRs' 
giant flares. 



^Schwartz (2005) measured the e-folding rising time to be ~ 4.9 ms from the observations on giant flare 
from SGR 1860—20 with Chinese Double Star polar spacecraft. 
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3.2. QPOs in the flare tail phase (t > 20 s) 

In the tail phase of the giant flare, declining of the radiation flux is explained as the 
shrinking of the flreball surface area (Thompson & Duncan 1995). After the flreball has 
evaporated to somewhat a smaller size, the plasma left from the flreball may form flux 
tubes above the flreball, and the slow sausage oscillations of such tubes may cause the 
QPO phenomena. The pulse phase-dependence of the QPO amplitude indicate that these 
oscillations are intrinsic to the neutron star surface (Israel et al. 2005; Strohmayer & Watts 
2006). This could be naturally explained as that the magnetic flux tube's footpoints are 
anchored at certain regions of the neutron star's crust. Sketche of the structure of such flux 
tubes is shown in Fig. [TJ The length of the tube L is related to its height H as 

L ^ hTTH (4) 

with the geometry index 6 ~ 1. Combine Eqs. (2) and Eqs. (4), we have the oscillation 
frequency in the slow standing sausage mode as 

where Hq = H/10^ cm, and the term in the bracket is attributed to the gravitational redshift 
effect (Mns and i?NS are neutron star mass and radius, respectively). 

We then choose typical values of T and H in the tail phase of the giant flare to estimate 
the QPO frequencies. In Feroci et al. (2001), a blackbody component with T = 9.3 keV, 
which accounts for ~ 85% of the total energy released above 25 keV, was derived at t = 
65 — 195 s after the onset of the SGR 1900+14 giant flare. So we may set kT ~ 10 keV for 
the flare in the pulsating tail phase (Thompson & Duncan 1995; Thompson & Duncan 2001; 
Hurley et al. 2005; Boggs et al. 2007). Previous studies of giant flares suggested a typical 
length scale L ~ 10 km for the 'flreball' formed in the burst (Thompson & Duncan 1995; 
Thompson & Duncan 2001; Hurley et al. 2005; Boggs et al. 2007), so we may take if g ~ 1 for 
the flux tube. These values give /slow ~ 19 Hz for the fundamental mode (j = 1) in a neutron 
star of mass Mns = 1.4 M© and radius R^s = 10 km. Consider its harmonic oscillations, 
most of the observed QPOs (summarized in Table 1) may be well explained, except those 
with very high frequencies of 625 Hz and 1840 Hz, because excitations of oscillations with 
j > 20 must be very difficult. 



3.3. Modulation of the Radiation 

In last section we showed that the slow sausage mode oscillations of flux tubes seem 
to be consistent with the QPOs observed. Now we move to the question how such tube 
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oscillations modulate the radiation observed in the tail phase of the giant flare. The mass of 
the plasma in the flux tube can be estimated to be AM ~ 10^^ g from Eq. (22) in Thompson 
and Duncan (1995), assume the total giant flare energy £" ~ 3 x 10^^ ergs (Cameron et al. 
2005). The plasma density is then p = AM/AV ~ lO^^/lO^^ ~ 10^ g cm^^ where AV is 
the tube volume. The Rosseland mean scattering cross-section in the direction parallel to 
the magnetic field is (Thompson & Duncan 1995) 

ffcs = 2.2 X IO^T^^-Vt ~ 1.5 X IQ-^^T^B-^, (6) 

where (Tt = (87r/3)(e^/'mcC^)^ is the Thomson scattering cross-section. Then the optical 

depth of the tube is 

r = na^sL = — a^sL ~ IO^T^q^^^B^^Lq, (7) 

flip 

suggesting that the flux tube is optically thick. Since the flreball itself is also optically thick 
(Thompson & Duncan 1995), the fraction of the radiation flux from tube in the total flux 
from the flreball can be simply expressed as the ratio of their surface area perpendicular to 

the line of sight 

-^tubc 'S'tubc 2aL 2a ,^, 



n^ 



-'^fireball '-'fireball ^^fircball ^fireball 

if we assume that the thermal temperatures of the flreball and the tube are roughly the same 
(about 10 keV). Here a is the radius of the flux tube, rgrcbaii is the radius of the flreball, and 
L ~ vrrfirebaii- As the flux tube oscillates, its cross-section and surface area vary, and so does 
the thermal emission from the tube. The amplitude of the QPOs can be derived to be 

SFtuhc Sd Stuhe Sa 2a , , 

r\^ rsj _ ( y ) 

-^fireball a 5'fircball ^ rgrcball 

This amplitude is consistent with the observational values ~ 10% — 20% (Israel et al. 2005; 
Strohmayer & Watts 2006; Watts & Strohmayer 2006), only when the radial flux tube 
oscillation amplitude 6a is in the range (5a/rfij.ebaii ~ 5% — 10%. This may also explain why 
the QPOs emerged from the light curve after about ~ 100 s from the onset of the giant flare: 
at the very beginning of the giant flare, the flreball was so big that (rfircbaii is very large), the 
tube radiation was too weak compared with that from the flreball, to produce a detectable 
QPO. 

It is noted that the QPOs are more likely to be detected in hard X-ray band (Strohmayer 
& Watts 2005; Watts & Strohmayer 2006), and their amplitudes in soft X-ray band are not 
as strong as in hard X-rays. One example is the 84 Hz QPO from SGR 1900+14 (Strohmayer 
&Watts 2005): Its amplitude increased from < 14% in the < 18 keV band, to (20 ± 3)% 
in the 12 — 90 keV band, and (26 ±4)% in the > 30 keV band. There are several possible 
reasons for this QPO amplitude-energy band dependence. The flrst is the photon splitting 
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mechanism (Thompson & Duncan 1995). Diffusion of pliotons from tlie fireball is primarily 
in E-mode due to the different scattering cross-sections between the two-polarization modes. 
Before the E-mode photons reach the flux tube, the photons with energy higher than 40 
keV still suffer serious photon splitting effect in the magnetic fields higher than the quantum 
magnetic field (the QED field -Bqed = 4.4 x 10^^ Gauss). This will produce an excess in 
the 10 — 20 kev in the spectrum from the fireball (Lyubarsky 2002), decreasing the ratio 
of the photon fluxes from the tube and from the fireball in this energy band and hence the 
QPO amplitude. The second is that the cross section of electron scattering decreases with 
decreasing photon energy, so the low-energy photons seen by us come from deeper region in 
the fireball, where temperature is higher (see Fig. 2 in Ulmer 1994). This may also increase 
the low energy photon flux from the fireball and reduce the amplitude of the QPOs in the 
same energy band. The third is cyclotron scattering of thermal photons in the flux tube, 
which occurs at energy (Ho & Lai 2001; Zane et al 2001) 

^ ZeB ,Z,, B ,^ 

where A and Z are the mass number and charge number of the ion, respectively. If we adopt 
B ~ 10^^ Gauss and A = Z = 1 ioi proton cyclotron resonance, this energy is ~ 10 keV. 
Such proton resonant scattering will up-scatter the photons in soft X-ray energy from the 
flux tube to higher energy (Lyutikov & Gavriil 2006), thus increase the QPO amplitude in 
the hard X-ray band. 



3.4. Excitation of the QPOs 

In this subsection we derive the energy needed to excite the oscillations of the flux tube, 
and the constraints on the possible energy source. From Wang et al. (2003), the energy 
needed to excite the slow sausage mode oscillations can can be estimated to be 

AE ~ -AMv^ ~ iAM(^-^)2 ~ AMkT/m^ ~ 10^^ ergs, (11) 

where AM ~ 10^^ g and kT = 10 keV are the mass and temperature of the plasma within 
the flux tube. 

According to recent investigations on the excitation of standing slow mode oscillations 
in a flux tube (Taroyan et al. 2004, 2005), the excitation energy should be injected into 
the tube within a timescale similar to the oscillation period r if the oscillations are excited 
by the energy deposition at the footpoint of the tube. If the energy deposition is through 
thermal conduction (Spitzer 1962), the deposited energy would be AE ~ F^St, with F^ ~ 



(1.84 X 10 ^T^/^VT)/ In A ergs ^cm ^ being the heating flux and S the heating area. Taking 
T ~ 10^ K, VT ~ 10^ Kcm-\ S ~ 10^ cm^, and r ~ 0.01 s, we get AE ~ 10^^ ergs, which is 
far less than that needed to excite the oscillations. So we conclude that the excitation energy 
is unlikely to be deposited through thermal conduction, but most likely by flare activities at 
the footpoint(s) of the tube. This may explain why there are no QPOs detected during the 
quiescent phase of SGRs. 



4. Discussion 

We have suggested the slow sausage oscillation modes of magnetic flux tubes to explain 
most of the QPOs detected in SGRs during giant flares. These QPOs generally last several 
rotational cycles, i.e. tens of seconds, and present useful constraints on the damping mech- 
anisms for the QPOs. The damping of magnetic loop oscillations in the solar corona has 
been extensively studied (Gaily 1986; Nakariakov et al. 1999; Ruderman & Roberts 2002; 
Stenuit et al. 1999; Taroyan et al. 2004). Radial wave leakage and resonant absorption are 
considered to be the main damping mechanisms. In our case, we consider only the effect of 
resonant absorption, as Ruderman & Roberts (2002) pointed out that, if the mass density 
inside the tube is great than that outside, the wave leakage effect is not important. Gonsid- 
ering the curvature of the tube, we can use the damping time scale of kink mode oscillations 
as that of sausage mode oscillation (Roberts 2000), which is given by (Ruderman & Roberts 

2002) 

2a po + p 

Td = — r. (12) 

7il po - Pe 

Here we assume that the mass density varies in the annuls region a — I < r < a from po to 

Pe- From Eqs. (18) and (22) in Thompson & Duncan (1995), we have a ~ 27^/ ~ 10^/. So 

Td ~ lO'^r, which is roughly consistent with the observational result r^ ~ 5 x lO^r for the 92 

Hz QPO in SGR 1806—20 lasting about 50 s. Note that the damping timescale is proportional 

to oscillation period r, so higher frequency oscillations should decay more quickly, which 

is also compatible with observations (Strohmayer 2007). However, the above results are 

based on the assumptions that the flux tube is thin, axis- symmetric with homogeneous mass 

distribution, which may not be satisfled in the real situation. Non-ideal effects, like the 

density stratiflcation, magnetic fleld curvature and twist, and thick tube limitation have 

been studied in both theories and numerical simulations (Van Doorsselaere et al. 2004a, b; 

Andries et al. 2005; Robert & Viktor 2006; Arregui et al. 2007). These works show that the 

above effects may only change the oscillation frequency by as much as 10% — 15%, but may 

seriously damp the oscillation mode and reduce the damping time scale. So the damping 

time scale with Eq. (12) should be taken as an upper limit 
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Observations indicate the QPOs frequencies may evolve with time. For example, Israel 
et al. (2005) found a possible time evolution of the QPO in SGR 1806—20 with the frequency 
increasing from 92.5 Hz to 95 Hz. From Eq. (2) it is seen that the oscillation frequency is 
determined by the length and temperature of the flux tube. Since the flux tube's footpoints 
are anchored at the surface of the star and confined by the super-strong magnetic fields 
that dominate the plasma's motion, its structure and length may not change much. The 
temperature of the flux tube is likely to remain nearly unchanged during the life time of the 
QPOs (less than tens of seconds), if there is no extra energy injected into the flux tube. So 
we would not expect a considerable change in the QPO frequency. The QPO amplitudes 
depend on the surface area of the fireball and the oscillation amplitude of the flux tube (see 
Eq. (9)), both of which decrease with time during the tail phase of the giant flare. The fact 
that the QPO's lifetime (< 50 s) is less than the evaporation time of the fireball (200 — 400 
s) implies that the tube oscillations damp faster than the shrinking of the fireball. So in this 
model the QPOs amplitude is predicted to decrease with time, which could be testified in 
the future high time resolution observations. 

In the tube oscillation model the oscillation frequency is /j oc j, while in the seismic 
vibration model /j oc [j{j + 1)]^^^, where j is the number of nodes (McDermott, van Horn & 
Hansen 1988 et al. 1988). To account for the observed frequencies with seismic vibration, 
one has to use the j = 2, 4, and 6 modes but ignore the j = 3 and 5 modes (Watts 
& Strohmayer 2006). Recently, Samuelsson & Andersson (2007) have derived a modified 
relation /j oc [(j — l)(j + 2)]^/^ for the torsional seismic modes using a general relativistic 
formulation, but the problem why the j = 3 and 5 modes do not exist still remains. In this 
point of view, the tube oscillation model seems to be more natural, since the j = 1, 2, and 
3 modes have all been used. 

There seem to exist three QPOs with different fundamental frequencies (18 Hz, 26 Hz 
and 30 Hz) in SGR 1806—20 (Watts & Strohmayer 2006). In the seismic vibration model, 
only one fundamental oscillation mode can exist (Israel 2005; Watts & Strohmayer 2006). 
However, in the tube oscillation model this could be explained if there are three flux tubes 
oscillating in slightly different fundamental modes. The existence of the multifrequencies 
might be due to the deviations in the temperature, density, and strength or topology of the 
magnetic field around the active regions in the magnetosphere. It is noted, however, that it 
is difficuh to explain the 625 Hz and 1840 Hz QPOs observed in SGR 1900+14 in the tube 
oscillation model. For these extremely high frequencies, the n = 1 and 3 torsional shear 
mode vibration (Piro 2005) might be more reasonable explanation. This also suggests that 
the QPOs in the giant flares may not be homogeneous, and may have different origins. 
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Tubel 




Tube 2 



Fig. 1. — Sketch of the flux tube formed in the tail of SGRs giant flare. The height of the 
tube H has also been indicated here. 
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SGR 1806-20 


SGR 1900+14 


numb 


er of nodes 


18, 26, 30 


28 




j = l 




53 




j=2 


92 


84 




j=3 


150 


155 




j=5 


625 




j^ 


=21(?) 


1840 




j^ 


=61(?) 



Table 1: QPO frequencies (in Hz) detected in the tails of the giant flares from SGR 1900+14 
and SGR 1806-20, together with the number of nodes j fitted in our flux tube oscillation 
model (see text for more information). 



